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ABSTRACT 

Aims. To compute the chemical evolution of the Galactic bulge in the context of an inside-out 
model for the formation of the Milky Way. The model contains updated stellar yields from mas- 



, sive stars. The main purpose of the paper is to compare the predictions of this model with new 

^""^ observations of chemical abundance ratios and metallicity distributions in order to put constraints 



on the formation and evolution of the bulge. 

Methods. We computed the evolution of several a-elements and Fe and performed several tests 



o 

i—{ ' by varying different parameters such as star formation efficiency, slope of the initial mass func- 

C/j ' tion and infall timescale. We also tested the effect of adopting a primary nitrogen contribution 

a ■ 

from massive stars. 



Results. The [a/Fe] abundance ratios in the Bulge are predicted to be supersolar for a very large 
range in [Fe/H], each element having a different slope. These predictions are in very good agree- 
ment with most recent accurate abundance determinations. We also find a good fit of the most 
recent Bulge stellar metallicity distributions. 

Conclusions. We conclude that the Bulge formed on a very short timescale (even though 
timescales much shorter than ~ 0. 1 Gyr are excluded) with a quite high star formation efficiency 
of V 20 Gyr"' and with an initial mass function more skewed toward high masses (i.e. x < 0.95) 
than the solar neighbourhood and rest of the disk. The results obtained here are more robust than 
previous ones since they are based on very accurate abundance measurements. 

Key words. Galaxy: bulge. Galaxy: evolution. Galaxy: abundances 

1. Introduction 

Galactic bulges are spheroidal systems that are found in the centre of most spiral galaxies, and usu- 
ally possess metallicity, photometric and kinematic properties that separate them from the disk 
components. Studying the evolution of the bulge of the Galaxy is of general interest because 
of the broad similarity of its integrated light to elliptical galaxies and other spiral bulges (e.g. 
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Whitford, 1978; Maraston et al., 2003). Bulges and ellipticals are located in the same regions of 
the fundamental plane (Jablonka et al., 1996); the Galactic bulge has a wide range in [Fe/H] (Rich, 
1988; McWilliam & Rich, 2004; ZoccaU et al., 2003; Fulbright et al., 2005) and enhancement of 
or-elements (McWilham & Rich, 1994; McWilliam & Rich, 2004; Rich & Origha, 2005); the age 
dispersion in the bulge at distances greater than 300 pc from the nucleus is observed to be narrow. 
Terndrup (1988) first argued for a globular cluster-like nature to the bulge population. Ortolani et 
al. (1995) compared the luminosity function of a metal rich globular cluster with that of the bulge 
field Baade's Window. Recently, Kuijken & Rich (2002) and Zoccali et al. (2003) showed that 
when the bulge field is decontaminated from disk foreground stars by proper motion cleaning or 
statistical subtraction, the remaining population is indistinguishable from an old metal rich glob- 
ular cluster. These studies form the basis for a growing consensus that the bulge is old and that 
its formation timescale was relatively short (< 1 Gyr). Very recently, medium- and high-resolution 
spectroscopy of bulge stars have been performed (Ramirez et al. 2000; Fulbright et al., 2006a, 
FMR06a), providing further support to a fast formation of the bulge. 

Finally, a short formation timescale for the bulge is also suggested on theoretical grounds by 
Ehnegreen (1999), who argued that the potential well of the Galactic bulge is too deep to allow 
self -regulation and that most of the gas must have been converted into stars within a few dynam- 
ical timescales. Moreover, Sarajedini & Jablonka (2005) suggest thai, since the differences in the 
metallicity distributions (MDs) of the Milky Way and M31 halos find no correspondence in those of 
their bulges, the bulk of the stars in the bulges of both galaxies must have been in place before any 
accretion event that might have occurred in the halo could have any influence on them, supporting 
a common scenario for the formation of bulges. 

Wyse (1998) also showed that the MD in our Galaxy is not consistent with a picture where 
the bulge is formed via accretion of satelhtes (see Searle & Ziim, 1978). There is now additional 
evidence that the bulge is not formed by sateUites similar to those observed at the present day in 
the halos of the Milky Way and M31. McWilham et al. (2003) find a systematically decreased 
abundance of Mn in the Galactic bulge, compared to stars in the Sgr dwarf spheroidal. So, even 
though Sgr does in principle reach high metallicity, its detailed chemistry is different. 

The baseline model for the chemical evolution of the bulge might well be the one-zone model 
(Searle & Zinn, 1978) with its instantaneous recycling and closed-box assumptions. Early numer- 
ical models (Arimoto & Yoshii, 1986) produced the observed wide abundance range. When su- 
pernova yields and detailed stellar lifetimes are incorporated, more extensive predictions become 
possible. This is the case of the Galactic bulge model by Matteucci & Brocato (1990, MB90) who 
predicted that the [a/Fe] ratio for some elements (O, Si and Mg) should be supersolar over almost 
the whole metaUicity range, in analogy with the halo stars, as a consequence of assuming a fast 
bulge evolution which involved rapid gas enrichment in Fe mainly by Type II SNe. At that time, 
no data were available for detailed chemical abundances; the predictions of MB90 were later con- 
firmed for a few a-elements (Mg, Ti) by the observations of McWiUiam & Rich (1994, MR94). A 
few years later, Matteucci et al. (1999, MRM99) modeled the behaviour of a larger set of abun- 
dance ratios, by means of a detailed chemical evolution model whose parameters were calibrated 
so that the metallicity distribution observed by MR94 could be fit. They predicted the evolution of 
several abundance ratios that were meant to be confirmed or disproved by subsequent observations, 
namely that a-elements should in general be overabundant with respect to Fe, but some (e.g. Si, 
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Ca) less than other (e.g. O, Mg), and that the [^^C/Fe] ratio should be solar at all metalUcities. They 
concluded that an IMF index flatter {x = 1.1-1. 35) than that of the solar neighbourhood is needed 
for the MR94 metaUicity distribution to be reproduced, and that an evolution much faster than that 
in the solar neighbourhood and faster than that of the halo (see also Renzini, 1993) is necessary as 
well. Whatever secular processes may be at work with respect to the formation and maintenance of 
the bar, the chemistry and stellar ages require that the stars formed early and self-enriched rapidly. 

Not all models of the bulge support these conclusions. Samland et al. (1997) developed a self- 
consistent chemo-dynamical model for the evolution of the Milky Way components starting from 
a rotating protogalactic gas cloud in virial equilibrium, that collapses owing to dissipative cloud- 
cloud colUsions. They found that self-regulation due to a bursting star formation and subsequent 
injection of energy from Type n supemovae led to the development of "contrary flows", i.e. alter- 
nate coUapse and outflow episodes in the bulge. This caused a prolonged star formation episode 
lasting over ~ 4 x lO' yr. They included steUar nucleosynthesis of O, N and Fe, but claim that gas 
outflows prevent any clear correlation between local star formation rate and chemical enrichment. 
With their model, they could reproduce the oxygen gradient of H n regions in the equatorial plane 
of the Galactic disk and the metaUicity distribution of K giants in the bulge (Rich, 1988), field 
stars in the halo and G dwarfs in the disk, but they did not make predictions about the evolution of 
abundance ratios in the bulge. 

Molla et al. (2000) proposed a multiphase model in the context of the dissipative collapse 
scenario of the Eggen et al. (1962) picture. They suppose that the bulge formation occurred in two 
main infall episodes, the first from the halo to the bulge, on a timescale th = 0.7 Gyr (longer than 
that proposed by MRM99), and the second from the bulge to a so-called core population in the 
very nuclear region of the Galaxy, on a timescale tb » th- The three zones (halo, bulge, core) 
interact via supernova winds and gas infall. They concluded that there is no need for accretion of 
external material to reproduce the main properties of bulges and that the analogy to elUpticals is 
not justified. Because of their rather long timescale for the bulge formation, these authors did not 
predict a noticeable diff"erence in the trend of the [or/Fe] ratios but rather suggested that they behave 
more likely as in the solar neighbourhood (contrary to the first indications of or-enhancement by 
MR94). 

A more recent model was proposed by Costa et al. (2005), in which the best fit to the observa- 
tions relative to planetary nebulae (PNe) is achieved by means of a double infall model. An initial 
fast (0.1 Gyr) collapse of primordial gas is followed by a supernova-driven mass loss and then by 
a second, slower (2 Gyr) infall episode, enriched by the material ejected by the bulge during the 
first coUapse. Costa et al. (2005) claim that the mass loss is necessary to reproduce the abundance 
distribution observed in PNe, and because the predicted abundances would otherwise be higher 
than observed. With their model, they are able to reproduce the trend of [O/Fe] abundance ratio 
observed by Pompeia et al. (2003) and the data of nitrogen versus oxygen abundance observed by 
Escudero & Costa (2001) and Escudero et al. (2004). It must be noted however that Pompeia et 
al. (2003) obtained abundances for "bulge-like" dwarf stars. This "bulge-like" population consists 
of old (~ 10-11 Gyr), metal-rich nearby dwarfs with kinematics and metaUicity suggesting an 
inner disk or bulge origin and a mechanism of radial migration, perhaps caused by the action of 
a Galactic bar. However, the birthplace of these stars is not certain, therefore we decided to omit 
these data from our model discussions, preferring to consider those stars for which membership in 
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the present day bulge is secure. Moreover, as we shall see, the use of nitrogen abundance from PNe 
is questionable, since N is known to be also synthesized by their progenitors and therefore it might 
not be the pristine one. 

In this paper we want to test the hypothesis of a quick dissipational collapse via the study of 
the evolution of the abundance ratios coupled with considerations on the metallicity distribution. 
Ferreras et al. (2003) already tried to fit the stellar metallicity distributions of K giants of Sadler 
et al. (1996), Ibata & Gilmore (1995) and Zoccali et al. (2003), which are pertinent to different 
bulge fields, by means of a model of star formation and chemical evolution. Their model assumes 
a Schmidt law similar to that of the disk, and simple recipes with a few parameters controlling 
infall and continuous outflow of gas. They explore a large range in parameter space and conclude 
that timescales longer than ~ 1 Gyr must be excluded at the 90% confidence level, regardless of 
which field is being considered. We want to show that abundance ratios can provide an independent 
constraint for the bulge formation scenario since they difl'er depending on the star formation history 
(Matteucci, 2000). 

We focus our attention on the evolution of the a-elements, carbon and nitrogen as a function 
of metalUcity. New data have lately become available for the abundance ratios of these nucUds: a- 
elements in particular are of paramount importance in probing the star formation timescale, since 
the signature of a very short burst of star formation must result into an a-enhancement with respect 
to iron, whereas ~solar abundances would imply that Type la Supernovae (SNe) had time to pollute 
the ISM with iron-rich ejecta (see e.g. Matteucci & Greggio, 1986). Moreover, we also analysed 
the evolution of nitrogen and studied the possibility of primary N being produced by all massive 
stars at any metallicity, as described in Matteucci (1986). This seems to provide a good fit to 
[N/O] and [N/Fe] in the solar vicinity (Chiappini et al., 2005; Ballero et al. 2006; see also Meynet 
& Maeder, 2002). The present paper does not primarily mean to point at a particular model as the 
"best model", but rather has the aim of showing which are the possible effects of varying the above- 
mentioned parameters; in fact, since the results for the MD of bulge giants are still preliminary, and 
high-resolution data are awaited, it is premature to draw firm conclusions on which combination 
of parameters provides the best fit. However, we can quite safely restrict the ensemble of plausible 
models by means of the present analysis. 

The paper is organized as follows: §2 describes the chemical evolution model, §3 provides 
a review of the data employed, §4 shows the model results regarding the supernova rates, the 
chemical evolution of abundance ratios, the MD of G, K and M giants in the bulge and in §5 we 
draw our conclusions. 

2. The chemical evolution models 

The adopted basic chemical evolution model closely follows that in MRM99. The main assumption 
is that the Galactic bulge formed with the fast collapse of primordial gas (the same gas out of 
which the halo was formed) accumulating in the centre of our Galaxy. We recall the fundamental 
ingredients of this model: 

- Instantaneous mixing approximation: the gas over the whole bulge is homogeneous and well 
mixed at any time. 
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Star formation rate (SFR) parametrized as follows: 

il,{r,t)^vG'{r,t) (1) 

where v is the star formation efficiency (i.e. the inverse of the timescale of star formation) in 
the bulge, A: = 1 is chosen to recover the star formation law employed in models of spheroids 
(e.g. by Matteucci, 1992) and G(r,t) - o-gas(r,t)/o-(r,tc) is the normalized gas surface mass 
density (where (Tgasir, t) is the gas surface mass density and cr(r, to) is the surface gas density 
of the bulge at the present time Iq - 13.7 Gyr). We tested also different values for k (such as 
1.5, which is the disk value) and the result do not differ much. The main difference with the 
solar neighbourhood resides in the higher v parameter for the bulge. 

We did not adopt a thresold surface gas density for the onset of star formation such as that 
proposed by Kennicutt (1998) for the solar neighbourhood, since it is derived for self-regulated 
disks and there seems to be no reason for it to hold also in early galaxy evolutionary conditions 
and in bulges. However, we also checked that adopting a threshold of 4 Mopc"^ such as that 
proposed by Elmegreen (1999) does not change our results, since a wind (see below) develops 
much before such a low gas density is reached. 

The initial mass function (IMF) is expressed as a power law with index x: 

(Pirn) oc m^^'^-^' (2) 

within the mass range 0.1 - lOOM©. 

The gas forming the bulge has a primordial chemical composition and accretes at a rate given 
by: 

GMnf^^^e-'l^ (3) 
o-(r, tc) 

where t is an appropriate collapse timescale and Air) is constrained by the requirement of 
reproducing the current total surface mass density in the Galactic bulge. Actually, we should 
use the halo chemical composition for the infalling gas, but it can be demonstrated that unless 
very high a-enhancements are adopted, the results are essentially the same. In principle, a 
slightly enriched infall could improve the fit of the MDs (see later) but this is not likely to have 
significant effects for realistic levels of enrichment. 

The instantaneous recycling approximation is relaxed; stellar lifetimes are taken into account 
in detail following the prescriptions of Kodama (1997). 

Detailed nucleosynthesis prescriptions are taken from Frangois et al. (2004), who made use of 
widely adopted stellar yields and compared the results obtained by including these yields in the 
model from Chiappini et al. (2003a) with the observational data, with the aim of constraining 
the stellar nucleosynthesis. Namely, for low- and intermediate-mass (0.8 - 8M0) stars, which 
produce '^C, N and heavy s-elements, yields are taken from the standard model of Van den 
Hoek & Groenewegen (1997) as a function of the initial stellar metallicity. Concerning massive 
stars (M > IOM0), in order to best fit the data in the solar neighbourhood, when adopting 
Woosley & Weaver (1995)Q yields, Frangois et al. (2004) found that O yields should be adopted 



' The models of Woosley & Weaver (1995) do not include mass loss, which might not be realistic for non- 
zero metallicity stars. However, it was shown (Bressan et al., 1981; Bertelli et al., 1990) that the combined 
effects of mass loss and overshooting compensate each other, in the sense that the resulting He-core mass 
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as a function of initial metallicity, Mg yields should be increased in stars with masses 1 1 -20Mq 
and decreased in stars larger than 20Mq, and that Si yields should be slightly increased in stars 
above AQMq, we use their constraints on the stellar nucleosynthesis to test whether the same 
prescriptions give good results for the Galactic bulge. Yields in the mass range 40 - lOOM© 
were not computed by Woosley and Weaver (1995), therefore one has to extrapolate them for 
chemical evolution purposes!^ 

- The Type la SN rate was computed according to Greggio & Renzini (1983) and Matteucci 
& Recchi (2001). Yields are taken from Iwamoto et al. (1999) which is an updated version 
of model W7 (single degenerate) from Nomoto et al. (1984). These supernovae are the main 
contributors of Fe and produce small amounts of light elements; they also contribute to some 
extent to the enrichment in Si and Ca. 

- Contrary to MRM99, we introduced the treatment of a supernova-driven galactic wind in anal- 
ogy with ellipticals (e.g. Matteucci, 1994). The bulge lies in the potential well of the Galactic 
disk and very massive dark halo that provides a high binding energy. Therefore, previous chem- 
ical evolution models supposed that a wind should not develop and moreover, the occurrence 
of winds did not seem to be suitable to reproduce the metallicity distribution of the Milky Way 
components (Tosi et al. 1998), especially for the disk. Moreover, Elmegreen (1999) sustained 
that the bulge potential well is too deep to allow for self-regulation and that the gas should be 
converted into stars in only a few dynamical timescales. However, this scenario was not tested 
quantitatively. We therefore supposed that the bulge is bathed in a dark matter halo of mass 
100 times greater than that of the bulge itself (i.e. Mdark = 2 x lO'^M©) and with an effective 
radius r^ark - lOOr^ = 200 kpc, where is the effective radius of the bulge (Sersic) mass 
distribution. To compute the gas binding energy Eb^gasif) we have followed Matteucci (1992) 
who adopted the formulation of Bertin et al. (1992). They analysed the properties of a family 
of self-consistent spherical two-component models of elliptical galaxies, where the luminous 
mass is embedded in massive and diffuse dark halos, and in this context they computed bind- 
ing energy of the gas. A more refined treatment of the Galactic potential well would take into 
account also the contribution of the disk; however, it is easy to show that the main contributors 
to the bulge potential well are the bulge itself and the dark matter halo. The condition for the 
onset of the galactic wind is: 

Eth,SN(tcw) - Ei, gas(tcw) (4) 

where Eth,sNit) is the thermal energy of the gas at the time t owing to the energy deposited by 
SNe (II and la) (see Matteucci 1992 and Bradamante et al. 1998). At the specific time tew (the 
time for the occurrence of a galactic wind), all the remaining gas is expeled from the bulge, and 
both star formation and gas infall cease. However, in all cases, the galactic wind occurs when 
most of the gas has already been converted into stars, and its effect on chemical evolution is 

(where all the heavy elements are produced) remains approximately the same as that in models without mass 
loss and overshooting. Additionally, nucleosynthesis with mass loss and rotation has so far been computed 

only for a few elements (e.g. Maeder et al., 2005). 

^ The authors are aware that the extrapolation process is problematic. However, the behaviour above this 
mass is not clear, since a supernova explosion may occur with a large amount of fallback. Moreover, it was 
shown (Franfois et al., 2004) that it is impossible to reproduce the observations at low metallicities in the solar 
neighbourhood if no contribution from stars in this mass range is considered. 
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negligible. Therefore, we shall limit ourselves to switching off star formation and infall when 
tew is reached. 

- We suppose that feedback from the central black hole is negligible. We shall test this hypothesis 
quantitatively in a forthcoming paper. 

As a fiducial model, we adopt the one with the following reference parameters: v - 20 Gyr 
collapse timescale r = 0.1 Gyr and a two-slope IMF with index x - 0.33 for M < IM© and 
X = 0.95 for M > IMq (Matteucci & Tornambe, 1987). The choice of such a flat IMF for the 
lowest-mass stars is motivated by the Zoccali et al. (2000) work, who measured the luminosity 
function of lower main-sequence bulge stars and derived the mass function, which was found to 
be consistent with a power-law of index 0.33 + 0.07. The IMF index for intermediate-mass and 
massive stars which is slightly flatter than that adopted by MRM99 in order to reproduce the MD 
of bulge stars from Zoccali et al. (2003, Z03) and Fulbright et al. (2006a, FMR06a; see §|3]and|43] 
for details) instead ot that from MR94. We explored a number of other possibilities by varying the 
model parameters in the following way: 

- Star formation efficiency: v from 2 to 200 Gyr"'; 

- For the IMF above IMq, we have considered the cases suggested by Zoccali et al. (2000, ZOO) 
in their §8.3, i.e. their case 1 (hereafter ZOO-1) with x - 0.33 in the whole range of masses, 
tiieir case 3 (ZOO-3) for which x = 1.35 for M > IMq (Salpeter, 1955) and their case 4 (ZOO-4) 
in which = 1 .35 for IMq < M < 2Mq and x = 1 .7 for M > 2Mo(Scalo, 1986). Our reference 
model corresponds to their case 2 with x - 0.95 for M > IMqQ therefore we call it ZOO-2. 
We recall that the fraction A of binary systems giving rise to Type la SNe is a function of the 
adopted IMF (see Matteucci & Greggio, 1986). Owing to the lack of information concerning 
the Type la SN rate in the bulge, we calibrate such a fraction in order to reproduce the Mannucci 
et al. (2005) estimate of the SN rate of an elliptical galaxy of the same mass. 

- Infall timescale: t from 0.01 to 0.7 Gyr; the latest hypothesis follows the suggestion by Molla 
et al. (2000) who assume a slower timescale for the formation of the bulge. We refer to the 
timescale th which they chose for the gas collapse from halo to bulge. 

Table [1] summarizes the features of the considered models^ 



3. The data 

We now turn to the datasets against which the models will be compared. These are re-normalized 
to the solar abundances of Grevesse & Sauval (1998) so that an artificial dispersion associated with 
the adoption of different solar abundance values is corrected for. 

' Actually, ZOO selected as "IMF 2" the one with x = 1 for M > IMq, but we performed calculations with 
X = 0.95, which is very similar, for comparison to the IMF with x = 0.95 of (Matteucci & Tornambe, 1987). 

* In the first column are indicated the names of the models, or their distinctive characteristic when no name 
was specified; in the second, third and fourth column are shown the IMF index, star formation efficiency and 
gas infall timescale respectively (* "Scalo" means x = 1.35 for 1 < M/Mq < 2 and x = 1.7 for M > IMq). 
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Model name/specification 


X (M > iMq) 


y (Gyr ) 


r (Gyr) 


Fiducial model (ZOO-2) 


0.95 


20.0 


0.1 


ZOO-1 


0.33 


20.0 


0.1 


ZOO-3 


1.35 


20.0 


0.1 


ZOO-4 


Scale* 


20.0 


0.1 


y = 2 Gyr-' 


0.95 


2.0 


0.1 


y = 200 Gyr' 


0.95 


200.0 


0.1 


T = 0.01 Gyr 


0.95 


20.0 


0.01 


T = 0.7 Gyr 


0.95 


20.0 


0.7 


SI 


0.33 


200.0 


0.01 


S2 


0.95 


200.0 


0.01 


S3 


0.33 


2.0 


0.7 


S4 


0.95 


2.0 


0.7 


S5 


0.33 


0.5 


1.25 



Table 1. Features of the examined models. 



Metallicity distribution 

Data for the [Fe/H] distribution of red giant branch (RGB) and asymptotic giant branch (AGB) stars 
in the bulge were taken from Z03, who provided photometric determination of metallicities for 503 
bulge stars. By combining near-infrared data from the 2MASS survey, from the SOFI imager at 
ESQ NTT and the NICMOS camera on HST, plus optical images taken with the WFI at ESO/MPG 
2.2m telescope within the EIS PRE-FLAMES survey, they constructed a disk-decontaminated 
(Mk,V - K) color-magnitude diagram (CMD) of the bulge stellar population with very large statis- 
tics and small photometric errors, which was compared with the analytical RGB templates in order 
to derive the MD. The advantage of this approach is that it allows determinations of metallicities 
of a great number of stars, although the relationship between the position in the CMD and the 
metallicity can be somewhat uncertain. 

Since the template RGBs are on the [M/H] scale (where M stands for the total metal abundance), 
in order to obtain the [Fe/H] distribution the a-enhancement contribution was subtracted in the 
following way (Zoccali, private communication): 

{ [M/H] - 0.14 if [M/H] > -0.86 
[Fe/H] = (5) 
[ [M/H] - 0.21 if [M/H] < -0.86 

This relation assumes that the a-elements in the bulge follow the abundance trends of globular 
clusters in the halo. They found that the resulting MD contained somewhat less metal-poor stars 
relative to a closed-box gas exhaustion model and that the G-dwarf problem (i.e. the deficit of 
metal-poor stars relative to a simple model) may affect the Galactic bulge even though less severely 
than in the solar neighbourhood (see e.g. Hou et al., 1998). 

In figure[T]we compare the photometric MD derived by Z03 to the spectroscopic one of MR94. 
The two distributions are broadly consistent, with a somewhat less pronounced supersolar [Fe/H] 
tail and a slightly sharper peak in the photometric case. However, the position of the high [Fe/H] 
cutoff for the photometric MD is wholly dependent on the metallicity assigned to the only two 
template clusters, namely NGC 6528 and NGC 6553 available in the high-metallicity domain. 
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Fig. 1. Comparison between the bulge MD derived by Z03 (solid histogram), FMR06a (dashed 
histogram) and that derived with the spectroscopic survey of MR94 (dotted histogram). 



Nevertheless, there are also indications (Zoccali, private communication) that the MD of the bulge 
stars, when achieved with high-resolution spectroscopy, is comparable to or narrower than that 
showed here, and that there are a very few supersolar-metallicity stars. 

In the same figure we also show the data from FMR06a, who performed a new high-resolution 
(R ~ 45000 - 60000) analysis of 27 K-giants in the Baade's Window (b = -4°) sample with 
the HIRES spectrograph on the Keck I telescope and, after determining their Fe abundances, they 
used them as reference stars in order to re-calibrate the K-giant data from Rich (1988). They found 
that the derived MDs are slightly stretched toward both the metal-poor and metal-rich tails with 
respect to those derived in previous works (Rich, 1988; MR94; Sadler et al., 1996; Zoccali et al., 
2003), although the overall consistency among these different MDs is reasonable. The metal-rich 
end probably extends to spuriously high abundances. 



a-element and carbon abundances 

- Abundances of O, Mg, Si, Ca, C and Fe for stars and clusters in the bulge are taken from Origlia 
et al. (2002), Origlia & Rich (2004), Origlia et al. (2005) and Rich & Origlia (2005). We refer 
to these abundances as the "IR spectroscopic database" hereafter. These datasets were obtained 
using the NIRSPEC spectrograph at Keck II, which allowed the determination of near-infrared, 
high-resolution {R ~ 25000), high signal-to-noise ratio (S/N > 40) echelle spectra and were the 
only ones to provide C abundances for bulge giants. They used the 1.6//m region of the spec- 
trum, corresponding to the H band. In all cases, abundance analysis was performed by means 
of full spectrum synthesis and equivalent width measurements of representative lines. Reliable 
oxygen abundances were derived from a number of OH lines; similarly, the C abundance was 
derived from CO molecular lines, whereas strong atomic lines were measured for Mg, Si, Ca, Ti 
and Fe. The data include observations of bright giants in the cores of the bulge globular clusters 
Liller 1 and NGC 6553 (Origlia et al., 2002, see also Melendez et al. 2003), Terzan 4 and Terzan 
5 (Origlia & Rich, 2004), NGC 6342 and NGC 6528 (OrigHa et al., 2005, see also Carretta et 
al. 2001; Zoccali et al. 2004) and measurements of abundances of M giants in Baade's window 
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(Rich & Origlia, 2005). The typical errors are of ±0.1 dex. The main considerations that were 
drawn from these abundance analyses are that a-enhancement is safely determined in old stars 
with [Fe/H] as high as solar, pointing toward early formation and rapid enrichment in both 
clusters and field, which are Ukely to share a common formation history. The [C/Fe] abundance 
ratios can be depleted up to a factor of with respect to the solar value, as expected because 
of the first dredge-up and possibly extra-mixing mechanisms due to cool bottom processing, 
which are at work during the evolution along the RGB, as also indicated by the very low (< 10) 
'^C/'^C abundance ratio (see also Origlia et al. 2003 and references therein). The analysis of M 
giants yielded abundances similar to those obtained with high-resolution optical spectroscopy 
of K giants; there is an apparent lack of supersolar-[Fe/H] stars, but the sample is too small to 
draw firm conclusions. 

- For O, Mg, Si and Ca we also included the abundance measurements of Fulbright et al. (2006b, 
FMR06b), who used the same specfra as in FMR06a, i.e. obtained the specfra of 27 bulge K 
giant stars at the Keck I telescope using the HIRES echeUe specfrograph with high resolution 
{R ~ 45000 - 67000) and high signal-to-noise ratio. The typical errorbars are of ~ 0.1 dex. 
The outcome of their analysis is that all elements produced from massive stars (i.e. a-elements, 
plus Na and Al) show enhancement in bulge stars relative to both Galctic thick and thin disk, 
although oxygen shows a sharply decreasing trend for supersolar [Fe/H], which is attributed 
to a metallicity-dependent modulation of the oxygen yield from massive stars. These results 
suggest that massive stars contributed more to the chemical enrichment of the bulge than to the 
disk, and consequently that the timescale for bulge formation were shorter than that of the disk, 
although they did not exclude other possibiUties (such as, e.g., an IMF skewed to high masses). 

- Finally, oxygen data from ZoccaU et al. (2006, Z06) were also taken. In this paper, Fe and O 
abundances for 50 K giants in four fields {b = -6; Baade's Window; Blanco b = -12; NGC 
6553) towards the Galactic bulge were derived; oxygen abundance was measured from the 
forbidden Une at 6300A. A high resolution (R = 45000) was achieved with FLAMES-UVES at 
the VLT. The typical errorbars are of ~ 0.1 dex. Also in this case, [O/Fe] is found to be higher 
in bulge stars than both in thick and thin disk, and supports a scenario where the bulge formed 
before the disk and more rapidly, with a formation history similar to that of old early-type 
galaxies. 

N vs. O 

Our nitrogen abundances in the bulge are derived from PNe. This may represent a problem, since 
while the measured oxygen abundance represents the true value of the ISM out of which the PN 
progenitor was formed, the observed N abundance has contributions both from the pre-existent 
nitrogen and from that produced by the star itself during its lifetime and dredged-up to the stellar 
surface. Moreover, these data come from emission fines which have a very complicated depen- 
dence upon several parameters (such as temperature, density and metal content) and assumptions 
(e.g. the photoionization model, from where the largest uncertainties come). Therefore it might 
be dangerous to employ these measurements to trace galactic chemical evolution. In principle, it 
would be possible to discriminate between enriched and primordial N by means of the C/N ratio 
(which is very different in the two cases, being much lower in the case of nitrogen-enriched stars). 
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Unfortunately, measurements for carbon in bulge PNe are available to date only for a very limited 
set of objects (Webster, 1984; Walton et al., 1993; Liu et al., 2001). Another possibility is to make 
use of symbiotic stars, i.e. presumably M giants in binary systems with a white dwarf or another 
hot companion. The envelope of the symbiotic star is photoionized from the hard UV radiation, 
leading to recombination line emission. Since the envelope has being observed, the abundances 
may be less evolved than in the PNe (Nussbaumer et al., 1988). 

It is possible that the N enrichment is not very dramatic especially in non-Type 1 PNe, which 
constitute about 80% of the PN population in the Galaxy (Peimbert & Serrano, 1980). Moreover, 
since the bulge presumably has not formed stars for a long time. Type I PNe (which have high- 
mass progenitors and are the most nitrogen-enriched) are not expected to be frequent. This is also 
confirmed by Cuisinier et al. (2000) who found quite low [N/O] ratios in their sample if compared 
to that resulting from self-enrichment. Moreover, Luna & Costa (2005) measured the N/O ratio of 
43 symbiotic stars towards the Galactic bulge and, as can be seen in their figure 5, the values of 
log(N/0) are consistent with those coming from studies of PNe. 

We employed the compilation of Gorny et al. (2004), who observed 44 PNe in the direction 
of the Galactic bulge with the aim of discovering Wolf-Rayet stars at their centre. The spectra 
were obtained with the 1.9-m telescope at the South African Astronomical Observatory, with an 
average resolution of 1000. Furthermore, Gorny et al. (2004) also merged their data with other 
published ones. Namely, they included the samples observed by Cuisinier et al. (2000), Escudero 
& Costa (2001) and Escudero et al. (2004). They obtained a total of 164 objects, among which a 
clear segregation of the subsamples is seen, due to the different selection criteria adopted to define 
each sample, and therefore none of them is truly representative of the bulge PN population. By 
merging the datasets, a more complete view of this population is achieved. Updated reddening 
corrections were applied to the objects from Escudero & Costa (2001) and Escudero et al. (2004). 
The merged sample was divided into two classes (according to the criteria listed in Stasinska et 
al, 1991), the first including those objects which are likely to be physically related to the Galactic 
bulge and the second containing the remaining objects which most probably belong to the disk. We 
only selected objects belonging to the bulge which had a clear detection of oxygen and nitrogen 
emission features; the resulting sample includes 103 objects. Errors in abundance derivations from 
both observational and theoretical uncertainties are typically 0.2 - 0.3 dex for [O/H] and can be 
even larger for [N/O]. 

In the future, when very high resolution IR spectroscopy (R * 100,000) becomes possible, N 
estimates could be also derived in giant stars from the faint CN Unes that are de-blended from 
stronger the CO and OH lines at high resolution. 

4. Model results 

4.1. Supernova rates 

Fig. |2] shows the predicted time evolution of the rate of Type II and Type la SNe in the Galactic 
bulge; the former, that die on small timescales, closely reflect the evolution of the SFR (for sim- 
phcity, we only show the case of different star formation efficiencies). The secondary peaks of the 
SNIa rates are software-related and are mostly due to the discontinuities in the adopted stellar life- 
times (Kodama, 1997) but do not affect the results concerning chemical abundances. The break in 
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Fig. 2. Evolution of the rates of Type la (left) and Type II (right) supernovae where different values 
of the star formation efficiency (middle panel). 



the Type II SN rates corresponds to the suppression of the star formation rate due to the achievent 
of the condition expressed in Eq. |2] which quite intuitively occurs earlier for flatter IMFs, higher 
v's and/or lower r's. However, even without a galactic wind, the Type II SN rate would become 
negligible at the same epoch, owing to the small amount of gas left in the bulge at that time. 

In those cases where the star formation is "bursty" (e.g. high star formation efficiency and quick 
formation timescale), the peak of the Type la SN rate (and therefore of Fe enrichment) can occur 
even before 1 Gyr, which is the timescale for Fe-enrichment in the solar neighbourhood. In fact, 
the time of occurrence of this peak is very sensitive to the underlying star formation history and 
can be as low as ~ 0.2 Gyr in the case of r = 0.01 Gyr or v = 200 Gyr"' (see also Matteucci & 
Recchi, 2001). Therefore, we expect that: 

1. The down-turn (i.e. change of slope) of the [a/Fe] vs. [Fe/H] plots, corresponding to the begin- 
ning of a sustantial Fe enrichment, will occur later, in general, at higher [Fe/H] values relative 
to the change in slope in the solar vicinity. This is due to the high star formation rate and the 
short timescale for bulge formation. This fact was pointed out for the first time by MB90. 

2. Where the IMF is flatter, the abundance trend of (T-elements vs. Fe is flatter as well, because 
SNela are fewer than in a Salpeter or steeper IMF. We also expect that in the ZOO-4 case, which 
favours the creation of SNIa progenitors, Fe production will be enhanced. 

3. For steeper IMFs (e.g. Scalo, 1986, or Salpeter, 1955), the nitrogen enrichment from low- and 
intermediate-mass stars is enhanced with respect to O, which mostly comes from massive stars. 
As a result, the [N/O] vs. [O/H] plot should lie above the one of our fiducial model. Moreover, 
N primary production from low- and intermediate- mass stars should be enhanced as well, 
and the flattening of the above mentioned plot towards the latest evolutionary stages should be 
prolonged. In general, any model that favours the contribution of massive rather than low- and 
intermediate- mass stars (i.e. high v, low t) should predict a downturn of the [N/O] vs. [O/H] 
plot if compared to our fiducial model, and vice versa. 
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4.2. Metallicity distribution of buige giants 

Fig. [3]shows the predicted MD for giant stars compared to the data from Z03 and FMR06c|fl 

The choice of a different IMF does not primarily influence the spread of the distribution but 
rather shifts its peak on the [Fe/H] scale: in general, flatter IMFs result in distributions peaked 
at higher values of [Fe/H]. Since the star-forming phase is relatively short, massive stars will be 
important in the Fe enrichment; if we adopt the Scalo (1986) exponent (model ZOO-4) rather than 
the Salpeter (1955) one (model ZOO-3) for stars more massive than 2Mq, the resulting MD will 
change dramaticaUy, moving the peak toward significantly lower [Fe/H]. Considerations on the MD 
allow us to exclude IMFs steeper than x = 0.95, whereas plots for flatter IMFs are all consistent 
with the observed distributions. In fact, the MD obtained with the very flat (x = 0.33) IMF of 
case ZOO-1 is almost undistinguishable from that calculated in case ZOO-2, showing that the MD 
progressively becomes almost insensitive to the flattening of the IMF under a certain value of the 
index. This is maybe suggestive of a "saturation effect" in the Fe enrichment from Type II SNe. 
Namely, a flattening of the IMF below a certain value of the index x does not produce a sensible 
increase of the number of massive stars, and moreover, the Fe mass ejected is assumed to be the 
same (i.e. does not increase as a function of mass) for all masses above 4QMq (see 

If we choose a moderate star formation efficiency (v = 2 Gyr"', comparable to that of the 
galactic halo), the predicted metallicity distribution is broadened and the number of both high- and 
low-metallicity stars is somewhat overestimated. The observed peak is poorly reproduced (both in 
height and position) for all of the MDs considered. In contrast, the extremely high value of v = 200 
Gyr ' yields an extremely narrow distribution, with an excessive peak height (reaching 0.68, but is 
truncated in the figure to preserve clarity). 

A change in the infall timescale affects mainly the spread of the distribution and only slightly 
the position of the peak. The low-metallicity wing is especially sensitive to the value of t. The 
model with t = 0.01 Gyr approaches a closed-box model, i.e. a model where all the gas is already 
present from the beginning (which would correspond to the limit t = 0). The gas soon reaches 
high densities and is consumed very rapidly. Thus, the number of metal-poor stars is overestimated 
and the predicted distribution extends below the observed low-metallicity tail. This confirms the 
considerations of Z03, i.e. an equivalent (although less important) manifestation of the G-dwarf 
problem occurs also in the Galactic bulge. We also calculated the MD resulting from the adoption 
of a slightly longer infall timescale, namely t - 0.7 Gyr, that corresponds to the timescale th for 
collapse from halo to bulge in the model of Molla et al. (2000). As the figure shows, the t = 0.7 
Gyr MD is at variance with observations, having a serious deficit of both metal-poor and metal- 
rich stars and too high a peak with respect to both measured distributions (again higher than 0.5 
but truncated in the figure). 

The main conclusions which can be drawn for variations of a single parameter are: 

1 . Changing the IMF slope has the general effect of shifting the peak of the MD towards lower 
metallicities for steeper IMFs, and vice versa. For x < 0.95 such variations become much less 
evident. We want to point out that continuous wind such as that of the Samland et al. (1997) 

^ We must remark that we did not make an attempt to convolve the predictions with uncertainties. Since 
the observed distributions are probably broader than the "true" distribution, due to random errors, and may be 
offset due to systematic errors, their comparison is likely to be affected. 
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Fig. 3. Predicted metallicity distributions in our bulge models for single parameter variations rela- 
tive to the reference model, with different choices of the IMF index (upper panel; see text for de- 
tails), star formation efficiency (middle panel), and gas infall timescale (lower panel). Our fiducial 
model is represented by the solid Une. Data are from Z03 (soUd histogram) and FMR06a (dashed 
histogram). 
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Fig. 4. Predicted MDs in models with v - 200 Gyr ' and r = 0.01 Gyr (left panel) and in models 
with V - 2 Gyr"' and t = 0.7 Gyr, with the adoption of two different IMFs (see text for details). 
Data are from Z03 (soUd histogram) and FMR06a (dashed histogram). 



model cannot have the same effects of flattening the IMF, since outflows lower the efl'ective 
yield, and thus if we had a continuous outflow we would need an even flatter IMF to reproduce 
the observed MD and abundance ratios. This was shown i.e. for the galactic disk by Tosi et al. 
(1999, flieir figure 7). 

2. Lower star formation efficiencies give rise to broader MDs and viceversa. The same is true of 
shorter infall timescales, which additionally increase the number of low-metallicity stars. 

3. The effect of both star formation efficiency and infall timescale on the position of the MD peak 
is negliglible. 

We now explore whether any other combinations of parameters can fit the observed MD. Since 
the position of the peak is essentially only determined by the IMF, combining IMFs other than those 
of models ZOO-1 and ZOO-2 with different values of the other parameters would not be useful to 
reproduce the required MD. On the contrary, the effects of infall timescale and star formation effi- 
ciency could compensate each other; therefore, we also considered the following "supplementary" 
(S) models: 

- Model SI: IMF like in model ZOO-1, v = 200 Gyx'\ t = 0.01 Gyr. 

- Model S2: same as in SI, but with IMF like in model ZOO-2. 

- Model S3: IMF Uke in model ZOO-1, v = 2 Gyr ', t = 0.7 Gyr. 

- Model S4: same as in S3, but with IMF like in model ZOO-2. 

In Fig. |4]are shown the resulting MDs for these supplementary models. 

It is clear that while combining a high value of v and a low value of t exacerbates the bulge "G- 
dwarf problem", a longer formation timescale can combine with a milder star formation efficiency 
to give rise to a MD compatible with observations, and this is what happened in some of the 
previous models (Samland et al., 1997; MoUa et al., 2000). However, such a degeneracy cannot 
be pushed above a certain range of values. We investigated longer timescales and found out that it 
becomes progressively harder to mantain the required MD both in height and position, since the 
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result of combining much longer timescales with efficiencies suitable to keep the required shape is 
to shift the distribution with the required shape toward lower metallicities, and vice versa, if we try 
to keep the predicted MD at the right position, its shape becomes too narrow, somewhat in contrast 
with the observed MDs. In any case, as we shall see, such a degeneracy is definitively broken when 
we take into account the evolution of abundance ratios (see next Section). 

For T > 2 Gyr, some star formation is predicted at detectable levels at the present time, contrar- 
ily to observations. 

Therefore, we conclude that good agreement with the observed MD is achieved only if the 
bulge formed on a short timescale, with a flat IMF and with a rather high star formation efficiency. 

4.3. Evolution of abundance ratios 
Oxygen 

Figure |5] shows the evolutionary plots of the [O/Fe] abundance ratio with metallicity, which is the 
most sensitive to the adopted parameters. The models are compared to the dara of the IR spectro- 
scopic database, FMR06b and Z06. Each panel iUustrates the effects of changing one of the model 
parameters as stated in §|2] The last panel also shows the results obtained with models S3 and S4. 
We remind that when the bulk of SNela begin to explode, the [O/Fe] ratio has a downturn and a 
knee in the [O/Fe] vs. [Fe/H] curve is expected. 

The fiducial model provides a satisfactory fit to the existing measurements and predicts the 
correct amount of oxygen enhancement, even though the slope of the predicted plot in the fiducial 
model is slightly flatter than the observations for supersolar [Fe/H] . However, as we shall see, the 
fit to this abundance plot cannot be improved without violating other constraints. The amount of a- 
enhancement crucially depends on the IMF index, and more in particular on the IMF of the Type II 
SN progenitors. The larger the number of massive stars, the higher the plateau at low metallicities. 
Model ZOO-1 overproduces oxygen, and is not consistent with data. Fig.|5]also allows us to exclude 
a Scalo (1986) exponent (model ZOO-4) for massive stars in the bulge, since the corresponding 
evolution of [c/Fe] vs. [Fe/H] lies well below the observed data points. A Salpeter (1955) IMF 
(model ZOO-3) cannot be excluded on the basis of these abundance ratios(; however, it was ruled 
out by the MD plot. 

Concerning the star formation efficiency, a value of v = 20 Gyr ' is consistent with the ob- 
served abundance ratios, whereas lower or higher values seem to be at variance with the (few) 
lowest-metallicity data. The model with v = 200 Gyr"' predicts values of [O/Fe] which are larger 
than those of the fiducial model at low metallicities, but then the gas is consumed very rapidly and 
star formation cease; this gives rise to a sharper decrease in the [O/Fe] ratio at higher metallici- 
ties. The predicted slope is compatible with observations, however the absolute amount of oxygen 
enhancement is slightly overestimated, and in any case this model was aheady excluded on the 
basis of the MD (Fig. [3). It is noteworthy that even though the star formation stops much earlier in 
time for higher v's, the model trajectory still spans the same range in [Fe/H] for all values of the 
star formation efficiency v. This is because at high values of the star formation efficiency, higher 
metallicities are attained in shorter times. 
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Fig. 5. Evolution of [O/Fe] vs. [Fe/H] in the bulge for different values of the IMF index (upper left 
panel; see text for details), star formation efficiency (upper right panel), and infall timescale (lower 
left panel). The lower right panel shows the results obtained combining a longer infall timescales 
and milder star formation efficiencies for two different IMFs (see text for details). Triangles repre- 
sent data taken from the IR spectroscopic database (field stars are symbolized by filled triangles, 
stars in globular clusters by open triangles. See |3]for specific references). The open and filled 
circles represent the low- and high-quality data, respectively, from Zoccali et al. (2006). Finally, 
the squares are the data from Fulbright et al. (2006b) and Zoccali et al. (2006). 



On the basis of the MD, we already saw that models with v = 2 and 200 Gyr ' must be 
excluded; considerations derived using the [O/Fe] abundance ratio provide a useful consistency 
check. We then suggest that a value of 20 Gyr ' can fit both the MD and the abundance ratios. 

Instead, changing the infaU timescale from 0.01 to 0.7 Gyr has almost no effect on the evolution 
of the [O/Fe] abundance ratio, and this holds also for other a-elements, if we exclude a small 
improvement of the agreement with data in the r = 0.7 Gyr case (which however, as well as the 
V = 200 Gyr ' case, is excluded by the MD plot). Therefore, there is no way of distinguishing 
between these models on the basis of considerations about abundance ratios, although we already 
excluded the cases of t « 0.1 Gyr and t > 0.1 Gyr which yielded MDs in contrast with the 
observed ones. 
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Fig. 6. Evolution of [a/Fe] vs. [Fe/H] in the bulge for O, Mg, Si, Ca for the fiducial model (solid 
line). We also plotted a solar neighbourhood fiducial line (dotted line) for comparison. Data for 
oxygen are the same as in Fig. |5] For the other elements, data are taken from Fulbright et al. 
(2006b, open squares) and, for Mg, from the IR spectroscopic database (field stars: filled triangles, 
clusters: open triangles. See text for detailed references). 

Finally, we discuss the results obtained with the "supplementary" models S3 and S4, whose 
MD was consistent with observations. We can see that it is necessary to assume a very flat IMF 
(model S3) to avoid underestimating the [O/Fe] ratio with the adoption of t = 0.7 and v - 2 Gyr 
and the predicted slope is much flatter than the observed one. For even longer timescales, flattening 
the IMF is no longer sufficient to fit the data. This is due to the fact that Type la SNe have more time 
to pollute the ISM with Fe causing the plot to turn down, falling below the observations. Therefore, 
formation timescales longer than ~ 1 Gyr must be excluded in any case. This is the consequence 
of the time-delay model for Type la SNe in the case of different star formation histories, as already 
pointed out by MB90 and Matteucci (2000). 

Other Qf-elements 

In figure|6]are shown the evolutionary plots of the abundance ratios of several a-elements (O, Mg, 
Si and Ca) to Fe versus Fe abundance, compared to the same datasets as [O/Fe], with the exception 
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of Z06. We also show a solar neighbourhood fiducial Une for comparison, calculated according to 
Frangois et al. (2004). A striking aspect of the predicted [a/Fe] vs. [Fe/H] relation is that the slope 
of the [a/Fe] ratios changes only at [Fe/H] ^ for the fiducial model, in agreement with the data 
and at variance with the solar vicinity where this occurs at [Fe/H] - -1.0. Indeed, it is evident that 
a star formation history and an IMF such as those suitable for the solar neighbourhood give results 
which do not agree with the bulge observations, the general trend being to severely underestimate 
the data especially at [Fe/H]> -1, when Type I SNe, whose formation is favored by the steeper 
IMF, start polluting the ISM considerably. 

It is worth noting that our fiducial model predicts that the slopes of [O/Fe] and [Mg/Fe] vs. 
[Fe/H] are different, the latter being less steep. This is mainly due to the fact that whereas mag- 
nesium is mainly produced by a restricted range of stellar masses (around 20 - 25Mq) oxygen 
is produced by a broad interval of stellar masses, (from 10 to lOOM©) the O/Fe production ra- 
tio increasing with stellar mass. This diff'erent trend is present also in the solar vicinity data (see 
Frangois et al., 2004) and is confirmed in the bulge by the observations of FMR06, which cover a 
larger range of metallicities than the IR spectroscopica database (i.e. [Fe/H] > and [Fe/H] < -1). 

Carbon and nitrogen 

Figure sQ and [8] show the evolutionary behaviours of C and N, which, in contrast to the a-elements 
are mainly produced by low- and intermediate-mass stars (and therefore on long timescales). The 
available dataset for carbon is limited to the measurements of the giant stars, which however are 
known to undergo mixing processes along the RGB, severly affecting the carbon abundances (in 
fact, we can immediately see that none of the models he beneath the observations, since carbon is 
likely to be depleted during the evolution along the RGB). Hence, the comparison between models 
and observations cannot allow to firmly constrain the former and to reach strong conclusions in 
general. The observations can only provide a lower limit to any plausible model. 

The behaviour of [C/Fe] at low metalUcities is largely affected by the assumed star formation 
efficiency, just as in the case of oxygen, whereas the curves do not show really discernible differ- 
ences above the solar metallicity. A high value of the star formation efficiency (v - 200 Gyr ') 
enhances the number of stellar generations and therefore accelerates the formation of intermediate- 
mass stars which are the main C producers, but again [C/Fe] falls below the predictions of the 
fiducial model due to fast gas consumption, as was the case for oxygen. The steeply falling trend 
of [C/Fe] until -solar [Fe/H] is due to the fact that the C/Fe production factor is decreasing with 
decreasing mass for high-mass stars. In this case too, abundance data for metal-poor bulge stars 
can confirm or refute the conclusions reached on the basis of the MD, even ffiough large values of 
V seem again excluded. 

A change in the IMF has a major impact on the contribution from low-mass stars. The plots 

show a bump at [Fe/H] ranging from 0.4 to 0.4, which corresponds to the time at which the first 

intermediate mass stars contribute to C enrichment, with the exception of the ZOO-1 IMF, for which 
the star formation ceases before ffie bulk of intermediate-mass stars have time to die. This bump 
is not visible in ffie solar neighbourhood, and its occurrence in the bulge is related to ffie strong 
metallicity-dependence of the adopted C yields for low- and intermediate-mass stars (Van den Hoek 
& Groenewegen, 1997), combined with the fact that in the bulge high metallicities are reached very 



20 



S.K. Ballero et al.: Bulge evolution 







1 1 1 1 1 1 1 1 1 














I 


- ZOO-1 

- ZOO-2 


ZOO 
ZOO 


-3 - 


; - 


■■■ y = 2 Gyr ' 

— y = 20 Gyr i 

- ■ y - 200 Gyr-i 


- 




0,5 












- 


:c/Fe] 















- 
















- 




-0.5 
-1 












— 
- 

- 




I 


— 1— 1 — \ — \ — — 1 — \ — \ — 1 — — 1 

T- 01 Gyr 




1 — h- 




1 1 1 1 1 1 1 i i i 
— ZOO-2 


— 1- 

- 






- T= 0.1 Gyr 








- S3 








T= 0.7 Gyr 






: — 


— S4 






0.5 












- 

































- 
















- 




-0.5 
-1 


— A 

1 1 1 1 1 1 1 1 1 1 


1 






1 1 1 1 1 1 1 1 1 1 








-1 








-1 








[Fe/H] 








[Fe/H] 





Fig. 7. Evolution of [C/Fe] vs. [Fe/H] in the Galactic bulge for various choices of the IMF (upper 
left panel) star formation timescale (upper right panel), and infall timescale (lower left panel). The 
lower right panel shows the results obtained combining a longer infall timescale and a milder star 
formation efficiency for two different IMFs (see text for details). Data sources are the same as 
in Fig. |5] Data are from the IR spectroscopic database (field stars: filled triangles, clusters: open 
triangles. See text for detailed references). 



early. The feature is more pronounced for larger values of x because low- and intermediate-mass 
star formation is favoured by these steeper IMFs. This is particularly evident from the difference 
between the plots adopting a Salpeter (1955) and a Scalo (1986) IMF, respectively. Before the 
onset of the intermediate-mass dominated regime, the amount of C enrichment increases for flatter 
IMFs due to the enhanced C production from massive stars; on the contrary, when the massive 
star dominated regime ends, C production is increased remarkably for steeper IMFs due to the 
enhanced contribution from low- and intermediate-mass stars. Again, the Scalo IMF looks rather 
implausible, since the predicted trend almost lies beneath the data points, and this is likely to mean 
that C is underproduced in this model. 

No appreciable change is seen when the adopted gas infall timescale is varied, with the excep- 
tion of a slightly different position of the occurrence of the bump. The plots resulting from models 
S3 and S4 are also shown. Bearing in mind the paucity of data below [Fe/H] ~ -0.5, the model 
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Fig. 8. Evolution of [N/O] vs. [O/H] in the Galactic bulge for various choices of the IMF (upper 
left panel) star formation timescale (upper right panel), and infall timescale (lower left panel). The 
lower right panel shows the results obtained combining a longer infall timescale and a milder star 
formation efficiency for two different IMFs (see text for details). The abundances measured from 
bulge PNe are taken from Gorny et al. (2004). 

S4 predicts a trend which passes through the cluster of data points at [Fe/H] ~ 0.0, but since as we 
already stated C is likely to be affected by stellar evolution, this model probably underpredicts the 
[C/Fe] ratio, again due to the Fe contribution from Type la SNe. Instead, model S3 leads to a larger 
C production and therefore keeps the consistency with the observations, the trend being only flatter 
than that of the fiducial model. 

Nitrogen is the only element for which our reference model cannot yield a satisfactory agree- 
ment with data: in fact, it tends to lie under the observational points, following their lower envelope. 
This is true also for the solar vicinity. As we can see from Fig. |8] there is no way of reproducing 
the average trend of the [N/O] abundance ratio with [O/H] for single-parameter variations unless 
we adopt parameters which have proven to lead to results at variance with other constraints: e.g. a 
very low star formation efficiency, which results into an underproduction of oxygen with respect 
to nitrogen and was dismissed on the basis of the resulting MD (and also to some extent on the 
basis of the [0/Fe] plot); or a steep IMF, favouring the formation of low- and intermediate-mass 
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stars (which are supposed to produce the bulk of N), which did not reproduce either the MD or the 
evolution of [a/Fe] vs. [Fe/H]. 

A good agreement is instead achieved with the S4 case, since it achieves the same effect of a 
low efficiency of star formation (i.e. oxygen underproduction) while the longer formation timescale 
does not influence the evolution of abundance ratios. Models S3 still somewhat improves the match 
with data relative to the fiducial model, but the flatter IMF favours O enrichment with respect to 
model S4. 

However, there is another way to obtain an acceptable fit to the observed abundances. In stan- 
dard chemical evolution models nitrogen production from massive stars is supposed to be purely 
secondary, i.e. N is created starting from seed nuclei of C already present in the gas out of which 
the stars were born. Also low and intermediate mass stars are supposed to produce N in a sec- 
ondary fashion but some primary N can be produced in intermediate mass stars during the third 
dredge-up episodes in conjunction with hot-bottom burning (Renzini & Voli, 1981; Van den Hoek 
& Groenewegen, 1997). The consequence of the secondary production is that the abundance of 
nitrogen should increase with metallicity in the earliest evolutionary phases, and that is what most 
evolutionary models of the Galaxy predict, not only for the bulge, but also for the solar neighbour- 
hood (Ballero et al., 2006; Chiappini et al., 2005). 

Therefore, we investigate what happens to the fiducial model if we assume that massive stars of 
all masses produce a constant amount (0.065Mo) of primary N at every metallicity. This hypothesis 
follows from the heuristic model of Matteucci (1986). The results are shown in Fig. |9]together with 
the plot of model S4. We see that this model seems to reproduce the average trend of the obser- 
vations better than the standard model and, although it adopts an ad hoc assumption, it is useful 
to understand that some mechanism of primary production of nitrogen is likely to occur at any 
metallicity in massive stars. For example, Meynet & Maeder (2002) calculated that stellar rotation 
can produce primary nitrogen in massive stars, and although Chiappini et al. (2003b) demonstrated 
that their rotation yields are insufficient to produce the observed trend in the solar neighbourhood, 
this hypothesis is the most promising one, as shown by Chiappini et al. (2006). We thus suggest 
that a continuous primary N production from massive stars is necessary at any epoch, in analogy 
with what is required in the solar vicinity (see Chiappini et al. 2005, for an extensive discussion on 
this point). 

From Fig.|9]it can be also noticed that model S4 differs from model ZOO-2 with primary N for 
[O/H] < - 1 , since while in the latter the [N/O] ratio is almost independent of [O/H] , in the S4 model 
the plot curves down sharply at low metallicities due to secondary production from massive stars, 
and at [0/H]~ -1 primary production from low- and intermediate-mass stars becomes dominant. 

5. Summary and conclusions 

In the present paper, we have computed chemical evolution models for the Galactic bulge where 
we have considered a scenario involving evolutionary timescales much faster than in the solar 
neighbourhood and in the halo, as the most probable one. 

We have adopted new chemical yields suggested by Francois et al. (2004) which provided the 
best fit to the solar neighbourhood abundance trends. 
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Fig. 9. Evolution of [N/O] vs. [O/H] (right panels) in the Galactic bulge in our fiducial model 
compared to a model where primary production of N from massive stars is assumed (Matteucci, 
1986) and to model S4. Data for N and O in bulge PNe are from Gomy et al. (2004). 

New observations, unavailable at the time of earUer chemical evolution models (e.g., Matteucci 
& Brocato, 1990; Renzini, 1993; Matteucci et al., 1999; MoUa et al., 2000) have been used to test 
these assumptions. The agreement between the new observations and our predictions is quite good, 
especially for the Fe abundance distribution and the [a/Fe] trends. 

Order-of magnitude changes in the main parameters determining chemical evolution (star for- 
mation efiiciency, stellar initial mass function, timescale of gas collapse) were applied in order to 
explore the consequences of such changes on the predicted stellar metalhcity distribution and trend 
of chemical abundances as a function of metallicity. 

Our main conclusions can be summarized as follows: 

- A short formation timescale combined with a high star formation efficiency fits both the ob- 
served metallicity distribution and chemical abundance ratios. This is typical of a star forma- 
tion history in a burst regime, i.e. strongly concentrated in the first stages of the lifetime of the 
system and vanishing very quickly, in analogy with elliptical galaxies. We suggest an efficiency 
of star formation of the order of 20 Gyr"'. However the assumption of a closed box must be 
given up since it gives rise to an excessive amount of low-metalUcity stars (even though not as 
seriously as in the case of the G-dwarf problem for the solar neighbourhood). A finite, though 
small, accretion time is instead required. We suggest that this timescale should be of the order 
of 0.1 Gyr. 

- There exists a sort of degeneracy between the gas infall timescale and the efficiency of star 
formation, in the sense that values of r longer than 0.1 Gyr can combine with values of v 
lower than 20 Gyr"' in order to match the observed metallicity distribution. However, such a 
degeneracy is broken when we consider also the evolution of [O/Fe] vs. [Fe/H] and in any case 
values of t longer than 1 Gyr make it impossible to fit the observed metallicity distribution as 
well. 
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- An IMF flatter than those suitable for the solar neighbourhood properties, such as that of Scalo 
(1998) or Weidner & Kroupa (2005), is necessary to reproduce the observational constraints. 
Namely, a value ofx = 0.95 for massive stars (even smaUer than that proposed by Matteucci et 
al., 1999, which was x = 1.1) gave the best overall fit. This can be theoretically understood if 
we note that the star formation in the Bulge proceeds like in a burst and there are several sugges- 
tions in the literature about a top-heavy IMF in starbursts (e.g. Baugh et al. 2005, Nagashima et 
al. 2005). Figer (2005) also finds a flat IMF in the Arches cluster near the Galactic center. The 
adopted flattening below IMg, as suggested by luminosity function measurements (Zoccali et 
al, 2000) does not aff"ect significantiy the abundance distribution. However, extremely flat IMFs 
(e.g. that of ZOO-1 model), if extrapolated to massive stars, wiU lead to some degree of oxygen 
overproduction. 

- The adoption of the fiducial model explains the behaviour of the difl'erent [a/Fe] abundance 
ratios with metalhcity very weU, and predicts difl'erent slopes for difl'erent or-elements, accord- 
ing to their nucleosynthesis. There is no need to invoke a second infall episode, as suggested by 
other authors (Molla et al., 2000; Costa et al., 2005) to explain the observed values of the abun- 
dance ratios. We do not exclude however a second infall episode on a much longer timescale 
such as that hypothesized by Molla et al. (2000) as it may help explain the presence of a very 
young stellar population confined to the very centre of our Galaxy, but it must involve a minor 
fraction of the bulge gas mass. 

- A certain amount of primary nitrogen from massive stars might be required in order to repro- 
duce the average trend of [N/O] vs. [O/H] with the fiducial model. The same conclusion was 
reached for the solar neighbourhood by Chiappini et al. (2005) and BaUero et al. (2006). The 
phenomenon of primary N production from massive stars as a result of rotation was studied by 
Meynet & Maeder (2002) and Maeder et al (2005) and although their yields are stiU not suf- 
ficient to explain the observed trend, this seems the most promising way (see Chiappini et al. 
2006). These considerations foUowed from observations of N and O abundances in planetary 
nebulae; in order to firmly assess this point, abundance measurements in stars which have not 
experienced nitrogen self-enrichment are required, especially at low metallicities. In alterna- 
tive, carbon measurements in the same planetary nebulae we have considered might allow an 
estimate of the amount of N synthesized by the PN progenitor. 

- In the near feature new sets of empirical metallicity and abundance pattern distributions in dif- 
ferent fields and based on high resolution spectroscopy will also become available. This will 
dramatically improve the chance of a more robust and quantitative comparison between the- 
ory and observations, with the ultimate goal of drawing the formation and chemical evolution 
history of the Galactic Bulge. 
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